A REMARK ON THE POTENTIALS OF OPTIMAL 
TRANSPORT MAPS 



PAUL W.Y. LEE 

Abstract. Optimal maps, solutions to the optimal transporta- 
tion problems, arc completely determined by the corresponding 
c-convex potential functions. In this paper, we give simple suffi- 
cient conditions for a smooth function to be c-convex when the 
cost is given by minimizing a Lagrangian action. 



1. Introduction 

The theory of optimal transportation starts from the problem of 
moving one mass to another in the most efficient way. Mathematically, 
the masses are given by two Borel probability measures /j, and v on a 
manifold M. The efficiency is measured by a cost function c : M x 
M — > K. and the problem is to find a Borel map which minimizes the 
following total cost: 



(1) / c(x,ip(x))dfi(x) 

J M 

among all Borel maps ip : M — > M which push [i forward to v. Here 
the push forward ip^/j, of a measure /x by a Borel map </? is a measure 
defined by ip*n{U) = n(^~ l (U)) for all Borel sets U CM. 

When the transportation cost c is given by minimizing a Lagrangian 
action, the existence and uniqueness of solution to the above problem 
is known. More precisely, assume that M is a compact connected 
manifold with no boundary and let L : TM — > K be a smooth function, 
called Lagrangian, which satisfies the following: 

• the second derivative ^2 is positive definite, 

• L is superlinear (ie. lim.|„|_ HX) ffl = 00 and | • | denotes the 
norm of a Riemannian metric). 

Let c be the cost function defined by 

(2) c(x,y)=M [ L( 7 (t),j(t))dt, 
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where the infimum is taken over all smooth curves 7(-) connecting x 
and y (ie. 7(0) = x and 7(1) = y). 

Assuming that the first measure /U is absolutely continuous with re- 
spect to the Lebesgue measure, Bernard-Buffoni (generalizing the ear- 
lier works of Brenier [6] in the Euclidean case and McCann [II] in the 
Riemannian case) proved the existence and uniqueness of solution to 
the above optimal transportation problem (see also [HI El E] for various 
extensions of this result and [I3J E] for a detail introduction to the 
theory of optimal transportation problem). 

In order to state the precise result, we need a few definitions. First 
let us consider the smooth function H : T*M — > R, called Hamiltonian, 
defined on the cotangent bundle T*M by the Legendre transform of the 
Lagrangian L: 

H(x,p) = sup \p(v) — L(x, v)]. 

veT x M 

Let H be the Hamiltonian vector field on T*M defined in local co- 
ordinates (xi, ...,x n ,pi, ...,p n ) by 

n ^ fdH dH \ 

and let 4>t '■ T*M — > T*M be the flow of the Hamiltonian vector field 
H. 

Another notion that we need is c-convexity of a function. A function 
/ : M — > R is c-convex if there is a function / : M — > R such that 

f(x) = sup [f(y) - c(x,y)}. 

It is known that any c-convex function is Lipschitz (in fact locally 
semi-convex) and hence differentiable almost everywhere. 

Theorem 1.1. j6j[lTJ|5] Assume that the measure \i is absolutely con- 
tinuous with respect to the Lebesgue measure and the manifold M is 
compact. Then the optimal transportation problem corresponding to 
the transportation cost (0|) has a unique solution ip : M — )■ M. More- 
over, there exists a c-convex function f : M — > R such that the solution 
ip is given by 

<p(x) = n((j)i(df x )), 
where n : T*M — > R is the canonical projection. 

The unique solution ip in Theorem ll.ll is called the optimal map. One 
natural question after Theorem 1 1 . 1 1 would be the following: Is there any 
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simple condition which guarantees a given function / : M — > R to be 
c-convex or whether the corresponding map <p(x) := 7r(<f>i(df x )) is an 
optimal map? 

In this paper, we give a solution to the above problem using curvature 
type invariants introduced in [2] (see also Section [2] and Theorem 13. II of 
this paper for more detail). We state a simple consequence of the main 
result (Theorem I3.ip when the Lagrangian is natural mechanical. More 
precisely, let (•, •) be a Riemannian metric on a compact manifold M 
and let U : M — > R be a smooth function, called potential. The next 
theorem is a version of Theorem 13.11 specialized to the Lagrangians of 
the form L(x, v) = ^\v\ 2 — U(x), called natural mechanical Lagrangians. 



Theorem 1.2. Assume that the Lagrangian L is of the form L(x, v) = 
t;\v\ 2 — U(x), the sectional curvature is non-positive, and the Hessian 
of the potential U : M — > R satisfies HessU < kl for some constant k. 
Let f be a C 2 function which satisfies 

(-y/\k\cath(y/\k\)I lfk<0 

Hessf > l-I ifk = 

[-y/\k\cot(y/\k\)I lfk>0. 

Then f is c-convex and the map <fi(x) := ir(<fii(df x )) is the unique opti- 
mal map pushing any Borel probability measure fi forward to ip*n. 



In the Riemannian case where the potential {7 = 0, Theorem 13. II can 
be improved using homogeneity of the corresponding Hamiltonian. Let 
,V2---,v n be an orthonormal basis of a tangent space T X M. 



|V/(as) 



We let S be the Hessian matrix of / with respect to this basis: 

n 

Ressf(vi) = ^S i: 



>ijVj. 



Theorem 1.3. Assume that the Lagrangian L is of the form L(x, v) = 
^\v\ 2 and the sectional curvature K < k for some constant k. Let f be 
a C 2 function and let X := 



|/c||V/|. Assume that f satisfies 




-Acoth(A)/ 



S > < -I 



-1 
-Acot(A)/ 



fc < 



k = 



k > 
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Then f is c-convex and the map tp(x) := 7i(<j)i(df x )) is the unique opti- 
mal map pushing any Borel probability measure \i forward to (p*/i. 

Note that if the manifold M has non-positive sectional curvature, 
then the condition in Theorem II .31 is just Hess / > — /. If the manifold 
is two dimensional and the Riemannian metric (•, •) is compatible with 
an almost complex structure J (ie. there is a endomorphism J : TM — > 
TM such that J 2 = —I and (Ju, Jv) = (u,v)), then the frame vi,v 2 
becomes t^4 and j^-, respectively. Therefore, Theorem 11.31 simplifies 
to the following. 

Theorem 1.4. Assume that the manifold M is two dimensional and 
there is a complex structure J which is compatible with the Riemannian 
metric (-,•)• Let L be the Lagrangian L(x,v) = ||v| 2 and assume that 
the Gauss curvature K < k for some constant k. Let f be aC 2 function 
and let X := y/\k\\V f(x)\. We define the following functions: 

{A coth(A) if k < 
1 ifk = 

Acot(A) ifk>0, 

h 2 = ((Hessf(JVf),JVf). 

Assume that the following inequalities hold for all points x for which 
Vf(x) + 0: 

. det Hessf > -| V/| 2 (^i + h 2 + £|V/| 2 ), 
. h 1 + h 2 >-(Z + l)\Vf\ 2 . 

Then f is c-convex and the map (p(x) := 7i(4>i(df x )) is the unique opti- 
mal map pushing any Borel probability measure fi forward to Lp*fi. 

The paper is organized as follows. In Section El various notions 
about curvature of Hamiltonian system needed in this paper is recalled. 
The general result (Theorem 13. ip mentioned above together with its 
consequence ( Theorem [L2]), are stated and proved in Section[3j Finally, 
the proof of Theorem 11.31 is given in Section HI 
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2. Curvature of Hamiltonian System 

In this section, we recall some definitions and properties about cur- 
vature of Hamiltonian system needed in this paper (see [21 [1] for a more 
detail discussion). 

Let H : T*M — > M. be a Hamiltonian and let H be the corresponding 
Hamiltonian vector field defined by 

u(H,-) = -dH(-), 

where u is the standard symplectic structure on the cotangent bundle 
T*M (ie. if (x±, ...,x n ,px, ...,p n ) is a local coordinate system of T*M, 
then co is given by u := J27=i ^Vi A dxi). 

Recall that a Lagrangian subspace of a 2n-dimensional symplectic 
vector space is a n- dimensional subspace on which the symplectic form 
vanishes. Let a be a point of the cotangent bundle T*M and let V a 
be the subspace, called the vertical space, of the tangent space T a T*M 
defined by 

(3) V a := {X e T a T*M\dn(X) = 0}, 

where 7r : T*M — >• M is the cotangent bundle projection. 

The vertical space V a is a Lagrangian subspace of the symplectic 
vector space T a T*M. If <j) t is the flow of the Hamiltonian vector field 
H, then we can define the following 1-parameter family of Lagrangian 
subspaces in T a T*M. 

J a (t) = dfc\V Ma) ) = {d<t>t\X)\X e V Ma) }. 

The family J a defines a curve, called the Jacobi curve at the point 
a, in the space of all Lagrangian subspaces contained in T a T*M. For 
each fix time t, we can define an inner product B l a on each subspace 
J aif) by 

B*(e,e) = -oj a (e,E(t )), 

where e(-) is a curve in the cotangent space T a T*M such that e(t) is 
contained in J a (t) for each t. 

The following proposition shows that if the Hamiltonian is fibrewise 
strictly convex, then the inner products B l are positive definite. 

Proposition 2.1. [T] Assume that the restriction H\t*m of the Hamil- 
tonian H to each fibre T*M is strictly convex. Then the bilinear form 
B l a is positive definite on J a (t). Moreover, B a is given by 
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For the rest of this paper, we assume that the Lagrangian L sat- 
isfies the conditions stated at the beginning of the introduction (ie. 
L is superlinear and the second derivative of L in the fibre direction 
is positive definite). Under these assumptions, the Hamiltonian H is 
fibrewise strictly convex and B l a are positive definite quadratic forms. 

Recall that a basis {e 1 , e n , f 1 , /"} in a symplectic vector space 
with symplectic form u is a Darboux basis if the following conditions 
are satisfied: 

u(e\ f) = 6a, u(e\ e>) = 0, u{f\ f) = 0. 

We can also pick a special family of basis e 1 (i), e n {t) of the subspace 
J a {t) which is orthonormal with respect to the inner product B a . More 
precisely, we have the following proposition (see [TJ for the proof). 

Proposition 2.2. [1] There exists a smooth family of basis e 1 (f),..., 
e n (t) on the vector space J a (t) orthonormal with respect to the canonical 
inner product £?* such that 

{e 1 ^),...,^),/ 1 ^) := -e\t),...,r(t) := -e n (t)} 

forms a Darboux basis of the symplectic vector space T a T*M. More- 
over, if (e 1 ^), e n {t)) is another such family, then there exists a con- 

n 

stant orthogonal matrix O such that e z (t) = ^ Oij&if)- 

i=i 

The 1-parameter family of Darboux basis 

{e 1 (t),...,e n (t)j\t),...J n (t)} 

in Proposition 12.21 is called a canonical frame at a. A canonical frame 
defines a splitting of the tangent bundle TT*M. More precisely, let 
% a be the subspace of the tangent space T a T*M, called the horizontal 
space at a, defined by 

(4) H a = spaQ{/ 1 (0),...,r(0)}. 

Then we have the splitting T a T*M = V a © H. a and both the vertical 
space V a and the horizontal a. are Lagrangian subspaces. 

Recall that n : T*M — > M is the canonical projection. The restriction 
of the differential dir to the horizontal space T-L a gives an identifica- 
tion between H a and T % ( a )M. We let v n be the unique vector in the 
horizontal space % a such that dir(v n ) = v. We call the vector v H the 
horizontal lift of v. 
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Let (x±, ...,x n ,Pi, ■■■iPn) be a local coordinate system of the cotan- 
gent bundle T*M and let v be a tangent vector of the manifold which is 
given by ^ Vid Xi in this system. Let Qj be structure constants defined 
by 

n / n 

v H = ^2v i ld Xi + ^2 Cijd P] 

i=l V J=l 

The following proposition gives a formula for the structure constants 

Cij. 

Proposition 2.3. [1] T/ie structure constants c^ satisfy 

2 Hp iPk CkiH PlPj 
k,l 

= y ] (-Hp k H PiPjXk — H Xk H PiPjPk — H PiXk H PkPj — H PiPk H XkPj ) . 

k 

The next proposition gives the definition of the curvature. 

Proposition 2.4. [I] Let {e 1 ^), e n (t), f n (t)} be a canoni- 

cal frame at a. Then there is a family of linear maps R a (t) : J a (i) — > 
J a {t) satisfying 

e\t) = -f\t) ) f = R a {t)e\t). 

The operators R a (t) are independent of the choice of canonical frame 
and it is symmetric with respect to the positive definite quadratic form 

B t a (R a (t)u,v) = B t a (u,R a (t)v). 

The operator R^ := R a (0) : V a — > V a is called the curvature operator 
of the Hamiltonian H. It also satisfies the following property 

(5) R a (t) = dfo'Rg^dfa. 

If X is a vector in the tangent space T a T*M, then the vertical part 
and the horizontal part of X are denoted by Xy and Xu, respectively. 
The following proposition gives a characterization of the curvature op- 
erator R^ without using canonical frame. 

Proposition 2.5. [T] Assume that X is a vertical vector field (ie. X(a) 
is contained in V a for each a ). Then the curvature operator R satisfies 

R%{V(a)) = -[H, [H,V] H } v (a). 
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Next, we consider the special case where the Hamiltonian is natural 
mechanical. More precisely, let (•, •) be a Riemannian metric on M. Let 
/ : TM — > T*M be the identification of the tangent bundle TM and 
the cotangent bundle T*M defined by I(v) = («,•). Let U : M E be 
smooth function, called potential. A Hamiltonian is natural mechanical 
if it is of the following form: 

H{a)= l -\I-\a)\ 2 + U{x), 

where x = ir(a) and it : T*M — > M is the canonical projection. 

The structure constants Cjj in this case is related to the Christoffel 
symbols as follows: 

Proposition 2.6. Assume that the Hamiltonian is natural mechanical 
and let be the Christoffel symbols corresponding to the given Rie- 
mannian metric and a local coordinate system (xi,...,x n ) of the man- 
ifold M. Then the structure constants c^j in the induced coordinate 
system (x\, x n ,pi, ...,p n ) of the cotangent bundle T*M are given by 

n 
k=l 

Proof. Assume that the Hamiltonian is given in local coordinates (x±, 
...,x re ,pi, ...,p n ) by 

H{x, p) = - ^2 9 l3 PiPj + U(x), 
hi 

where = \d Xi ,d Xj } and g l i denotes the inverse matrix of g^. 
By Proposition 12.31 the structure constants Qj is given by 

2 e = e U ,d -fp> - e-'T* - <> ik TA 

Tt \ dx k ox k dx k J 

s,r k,l x ' 

If we rewrite the above equation, then we get the following as claimed. 

□ 
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For each point a in the cotangent bundle, we can identify the cotan- 
gent space T*M with the vertical space V a by 

(6) a (-)■ a v := — (a + ta). 

at t=o 

We call a v the vertical lift of a at the point a. We also recall that 
/ : TM — y T*M is the identification of the tangent and the cotangent 
bundle induced by the Riemannian metric. Under these identifications, 
the curvature R H defined above and the Riemannian curvature 91 are 
related as follows. 



Proposition 2.7. Assume that the Hamiltonian H is natural mechan- 
ical. Then the curvature R H of the Hamiltonian H is given by 

Rf u ((Iv) v ) = (I(K(u,v)u + HessU(v))) v , 

where 9\ is the Riemannian curvature and Hess denotes the Hessian 
with respect to the given Riemannian metric. 



Proof. Let E : T*M — )■ R be the kinetic energy Hamiltonian defined 
by 

E{a)= l -\I-\a)\\ 
Since H = E + U, we have 

H = E - {dU) v . 

Let V be a vertical vector field. Since [(dU) v , V] is vertical, it follows 
from Proposition 12.51 that 

R H (V) = -[E-(dU) v ,[E,V] H ] v . 

By Proposition 12.61 the horizontal spaces of the Hamiltonians H and 
E coincide. It follows that 

R H (V) = R E (V) + [(dU) v , [E, V]n]v- 
By Theorem 5.1], 

RUilvY) = (m(u,v)u)) v . 
Therefore, it remains to show that 

[(dU) v , [E, {IX) V ] H ] V = (J(Hess[/pO)) v 
A calculation using local coordinates shows that 

dn({E,(IX) v ]) = -X. 
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Therefore, 

i(dU) V , [E, (IX) v } n } v = -[(dU) v , 1% = (/Hess U{X)) V . 

□ 



3. The Main Result 

In this section, we will state and prove the main result (Theorem 
13. ip . Before that, we need a few definitions. 

Let / : M — > R be a C 2 function. Its differential defines a map 
P '. x i — y df x from the manifold M to the cotangent bundle T*M. There- 
fore, the differential dP : TM — > TT*M of this map sends each tangent 
space T X M to a n dimensional Lagrangian subspace dP(T x M) of the 
tangent space T dfx T*M = V dfx © H dfx (see fl3} and (0} for the defini- 
tion of the vertical space V a and the horizontal space T-L a , respectively). 
Therefore, dP(T x M) is the graph of a linear map Hess 11 f : %df x Vdf x - 
We call this map the //-Hessian of /. When the Hamiltonian is nat- 
ural mechanical, the iZ-Hessian coincides with the usual Hessian de- 
fined using the given Riemannian metric (see Proposition 13. 3p . Let 
e 1 (t), e n (t), / 1 (t), f n {t) be a canonical frame at the point a (see 
Proposition 12.21 for the definition) and let S be the matrix representa- 
tion of the iZ-Hessian defined by 

n 

i=i 

Finally we can state the main result. 



Theorem 3.1. Assume that the curvature R H corresponding to the 
Hamiltonian H satisfies < kl for some constant k and all points a 
in the cotangent bundle T*M. Assume that f is a C 2 function which 
satisfies 

{-^\k\cot\l(yJ\k\)I lfk<0 

S > < -I if k = 

[-• v /pfc[cat(- N /pfe[).T ifk>0. 

Then f is c-convex and the C 1 map <p(x) := 7r(<pi(df x )) is the unique 
optimal map pushing any Borel probability measure /j, forward to tp*fi. 



Before giving the proof of Theorem 13.11 we note that Theorem 11.21 
is an immediate consequence of Proposition \2.7\ Theorem I3.1[ and the 
following two results (Proposition 13.21 and Proposition 13. 31) . 
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Proposition 3.2. Assume that the Hamiltonian is natural mechani- 
cal. Let Vi,...,v n be an orthonormal frame of the tangent space T X M. 
Then there is a unique canonical frame {e 1 (t), e n (£), f n {t)} 
satisfying 

e i (0) = (lv i ) v , f(0) = vl i , i = l,...,n. 



Proof of Proposition \3.2l Since Vi,...,v n is orthonormal, we have, by 
Proposition 12.11 

& a {(Iv i ) v ,(Iv j ) v ) = 6 ij . 

It follows from Proposition 12.21 that there exists a unique canonical 
frame {e 1 ^), e n (t), f\t), f n (t)} such that e*(0) = (Jw i ) v . 

Let (xi, x n ,pi, ...,p n ) be a local coordinate system of the cotan- 
gent bundle T*M around the point a. Assume that the vectors v 
and w are given in this coordinate system by v = Y^=i v i®xi and 
w = Y^i=\ w Ai, respectively. Then 

n n / n 

(Iv) v = ^2 QijVidxj, w n = Wj I d Xi + Cijd Pj | 

i=l i=l \ j=l 

where = (v^Vj). 
It follows that 

n 

u(v v , w H ) = 22 9ij v iWj = (v, w) . 

M=l 

Therefore, by the definition of e*(0) and the fact that v±, v n is an 
orthonormal family, we have 

u(e i (0),v?)=oj((Iv i j v ,v?)=5 ij . 

Finally since e 1 (0), e n (0), / 1 (0), f n (0) is a Darboux basis, we 
have v^- = / J '(0) as claimed. □ 



Proposition 3.3. Assume that the Hamiltonian is natural mechanical. 
Then 

Hess H f{v n ) = {I{Hessf{v))) v . 
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Proof of Proposition Iff. 31 Let (xi , . . . , x n , p\ , . . . , p n ) be a local coordinate 
system and recall that denotes the structure constants. Suppose v 
is given, in this coordinate system, by v — ^Vid Xi . By Proposition 
]6j its horizontal lift is given by 



5> (^+£r* Pfc *u 

i \ k,j / 



where denotes the Christoffel symbols. 

In the above coordinate system, the differential df is given by ^xi, 

Xn ' ~§ij m ^ f° nows from the definition of iFHessian that 

HeBsW) = f^T " E r ^) ^ = (/(Hess/(,))) v 

i,j V 1 J A; / 

□ 

Proof of Theorem \3.1[ Recall that (fit denotes the flow of the Hamil- 
tonian vector field H. We define the maps $t : M — > T*M and 

(ft : M M by 

$ t (x) = (j) t (df x ), <p t = Tro& t . 

The image At of the map c/$t : T^M — > T^ t ^T* M defined a La- 
grangian subspace of the symplectic vector space T$ t ( x )T* M . 
Let 

e\t),...,e n (t),f\t),...J n (t) 

be a canonical frame at the point &o(x). Let a t and fet be two family 
of matrices defines by 

n 

(7) W(0))) = + #e>"(f)). 

If we let Ft = (fi(t),...J n (t)) T and E t = ( ei (t), e„(t)) T , then © 
implies that at Ft + frf^t is constant in time t. By Proposition 12.41 if we 
differentiate the above equation with respect to t, then we get 

= a t F t + a t F t + b t E t + b t E t = a t F t + a t R t E t + 6^ - M 7 *, 

where R t denotes the matrix given by 

n 

(8) R t (e\t)) = Y,(Rt)^(t). 

3=1 
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It follows that 

(9) a t = b t , b t = -a t R t . 

Note that the Lagrangian subspace A , and hence A t for all small 
enough t, is transversal to the vertical space V$ t ( x ). This, in turn, is 
equivalent to the matrix a t being invertible for all such t. We claim 
that A t is transversal for all time t in the interval [0, 1]. In other words, 
we need to show that a t is invertible for all t in [0,1]. If we let 

(10) S t := a^b t , 

then it is enough to proof that S t is bounded for all time t in [0,1]. 
From (j7|), we have that 

n 

f (0) + Hess"/(n0)) = £(aJ"/»"(0) + #e»"(0)). 

It follows that a = I and Uess H f (f (0)) = YT j= i &o e*(0). Therefore, 
So satisfies the initial condition 

So = a 1 bQ = S. 

By (Q and (jlOp . the matrix St satisfies the following Riccati equation 

S t + S 2 t +R t = 0. 

Note that the matrix St is bounded above. For the lower bound, we 
need the following comparison principle of matrix Riccati equations. 
We denote the transpose of a matrix B by B T . 

Theorem 3.4. [TJ] Let S\ be the solutions of the matrix Riccati equa- 
tions 

si = Ai+Bisi+si(Bi) T +sicisi, 2 = 1,2. 

Assume that 

A^ B^ \ ( A 2 B^ \ , „i ^2 
i\t ru l — l / n2\T r*2 ) ana ^0 < ^o- 



Then S\ < Si 



Therefore, by the assumption and Theorem 13. 4[ if we consider the 
equation 

(11) B t + SI + kI = 0, 
then we have 

(12) S t < S t . 
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The solution to (TTTj) is given by the following theorem. 



Theorem 3.5. [10] Let S t be the solution of the matrix Riccati equation 
with constant coefficients 



St = A + BSt + S t D + StCSt. 



Let M t := ^ ^ 3 ^ 4 J be the fundamental solution of the following 
equation with initial condition M = I : 



B A 
-C -D 



z. 



Then S t = (MxS + M 2 )(M 3 S + M 4 ) 



It follows that the solution to the equation (jTTj) with initial condition 
So = S is given by 

s^r^ciyt))- 1 , 



where 



ri(t) 



cosh(y/\k\t)S 
S 

cos( J\k\t)S - 



k\sinh(y/\k\t)I iik<0 
if k = 

ia(J\k\t)I if k > 



and 



' sinh(W\k\t) 



r 2 (t) 



ts + i 

sin(-v/|fe|t) 



«S + cosh(-^|A;|t)7 if jfc < 
if fc = 

S + cos{J\k\t)I iffc>0. 



Therefore, by (jlSjl . & is bounded for all t in [0, 1] if r 2 (t) > 0. This, 
in turn, follows from the following 



(13) 



S > 



\k\coth(^\k\)I iik<0 
-I if k = 

rfcicot(vW) 7 iffe>o. 



This finishes the proof of the claim that Ft is transversal for all t 
in the interval [0, 1]. It follows from the claim and compactness of the 
manifold M that the map ip t is a diffeomorphism for each t in [0,1]. 

The following theorem is proved by the method of characteristics 
(see (U Theorem 17.1, Section 17.2] for a proof). 
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Theorem 3.6. Assume that tp t is a diffeomorphism for each time t in 
the interval [0,1]. Then the curve 

t ip t (x) 

is a strict minimizer of the minimization problem in (TJ|) for each point 
x. 

Moreover, there exists C 2 solution f t to the Hamilton- Jacobi equation 
d t f t + H(x,df t ) = 0, f = f 

and it satisfies 

3> t (x) = {df t ) Mx) . 

Finally, we show that ipi is the optimal map between any measure /i 
and (</?i)*/i. Let 7(t) be a minimizer of (EJ) which starts from 7(0) = x 
and ends at 7(1) = y. Then we have 

h(y) - f (x) = J j t ftd(t))dt = J fMt)) + df t (W))dt. 

By the Hamilton- Jacobi equation in Theorem l3.6[ the above equation 
becomes 

(14) h{y) - f (x) = [ -H(j(t), (d/ t ) 7W ) + df t (j(t))dt. 

Jo 

By the definition of the Hamiltonian H, the above equation gives 

(15) fi{y)-h{x)<l L(<y{t),i{t))dt = c(x,y). 



By Theorem 13 .6[ we have (ft(x) = d p H (tp t {x)). Therefore, if we 

dft 

let j(t) = (ft( x ), then fTl^ becomes 

fi(ipi(x)) - f (x) = / -H(x,(df t ) Vt{x) ) + dft{d p H ((p t (x)))dt. 

JO dft 



Finally, by the definition of the Hamiltonian H and Theorem 13.61 
the above equation gives 

fi((fi(x)) - f (x) = / L(<p t (x),<pt(x))dt = c(x,(p 1 (x)). 



If we integrate both sides with respect to fi, then we have 

(16) / c(x,<pi(x))d(ji= / fid(<p*(i)- / fodfi. 

Jm Jm Jm 
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Therefore, this finishes the proof if we combine ( fl5|) . ( fl6l) . and the 
following standard theorem in the theory of optimal transportation (see 
[H] for a proof). 

Theorem 3.7. Let 7Ti, 7^2 '■ M X M — > M be the projections onto the 
first and the second entry, respectively. Then the following holds: 

inf / c{x, y)dll = sup / fadv - / f dfi, 

JMxM JM JM 

where the infimum on the left is taken over all measures U on the 
product space MxM satisfying (ni)^U = /i and (7r 2 )*n = fi, and the 
supremum on the right is taken over all pairs of continuous functions 
(/o, fi) satisfying fi(y) - f (x) < c(x, y). 



□ 



4. The Riemannian case 

In the section, we specialize to the Riemannian case and give a proof 
of Theorem 11.31 In this case, we can use the homogeneity of the corre- 
sponding Hamiltonian to improve the result in Theorem 13.11 

Proposition 4.1. [21 Lemma 5.1] Assume that the Hamiltonian H is 
homogeneous of degree 5 + 1 in the fibre variable. Let f be the Reeb 
field defined by f(a) = a v . Then r(a) — t6H(a) is contained in J a (t) 
for all t. In particular, the Hamiltonian vector field H is horizontal. 



Proof of Theorem \1.3[ Let tp and $ as in the proof of Theorem 13. 11 Let 
e 1 (i), e n (t), f n {t) be a canonical frame at the point a = df x 

of the cotangent bundle T*M. We claim that e 1 (t) can be chosen to 

be z(t) := (r(a) - i#(a)) if a ^ 0. 

First, note that z(t) has norm one with respect to the inner product 
B f a . Indeed, by the definition of B f a , we have 

, / x /xx / / x ./xx u(r(a).H(a)) 

BMt),z(t)) = -U{z{t),z{t)) = |)_! Q) |2 ■ 



A calculation using local coordinates, we have 
u(f(a),H(a)) = I^V)! 2 
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and 

Bi(z(t),z(t)) = l. 

By Theorem 14.11 

z(t) = A t E h 

where A t is a n x 1 matrix and E t = (e (t), e n (t)) . 

If we differentiate the above equation twice and note that z(t) = 0, 
then we have 

A t E t - 2A t F t - A t R t E t = 0, 

where F t = (f(t), f n {t)f and R t is defined as in ©. 

It follows that At = and A t is independent of t. Therefore, by 
Theorem 12 A\ we can choose e 1 (t) to be z(t) by applying an appropriate 
orthogonal transformation if a = df x ^ 0. 

For the rest of the proof, we assume that e 1 (t) = ij-i( a )i (^( a ) ~ 

tH(a)) and = — e 1 (t) = |j-i( a )| -fl"(oO whenever 

Let be the components of the matrix R t defined as in Theorem 
13. 11 It follows that from the definition of the inner product and flSJ) 
that 

Rf = ^(^(e^)),^)) = -u a {d^ l R H Ma) dUe l {t)),r{t)). 

Since the symplectic form u is preserved along the Hamiltonian flow 
<f) t , we also have 

if = -u Ma) (R$ t{a) dMeKt)),dMf j m 

= Bl {a) {R» {a) dUe\t)),dUe j {*)))■ 

Let v l (t) be the tangent vectors defined along the geodesic t i— > tt o 
0t by the vertical lift: d<f> t (e l (t)) = (Iv l (t)) v . Note that since the 
frame e 1 ^), ...,e n (t), fk), f n (t) is a Darboux frame, v 1 ^), t> n (t) 
is orthonormal with respect to the Riemannian metric. It follows from 
Proposition 12.11 and 12.71 that 

^ = ^ (a) (i^ (a) (/t; i (t)) v .(M*)) V ) 

(17) =^ t(a) ((/0l(/- 1 (^(a)),i; i (t))J- 1 (^(a))) v ,(J^(t)) v ) 

= (5H(/- 1 (0 t (a)),^(t))/- 1 (^(a)),^(*)> 7r ^ (a)) • 

If we assume that V/(x) 7^ 0, then e 1 (t) = rFTT^yi (^(q) — tH(a)) 
and we have, by Proposition I4.1[ 

(18) (Iv\t)) v = dUe\t)) = - J—f^a)) = T7 ^ ^(«) V 
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Note that I 1 a = I l df = V/ and the Riemannian exponential map 
exp satisfies ir((p t (Iv)) = exp(w). Therefore, it follows from (fT7|) and 
(HID that 

(19) Ft} = \Vf(x)\ 2 (m(v\t),v\t))v\t),v\t)) cMtVf(x)) . 

Note that the equation in ffl9l) holds also in the case V/(x) = 0. 

By the assumption of the theorem, the sectional curvature is bounded 
above by k. Therefore, the following 

(9V(*),V(*),-) 

defines a bilinear form on the the orthogonal complement of v 1 (t) which 
is bounded above by kl. It follows from this and (TT9]) that 



(20) 



Rt < 











k\Vf\ 2 I 
Let S be the Hessian matrix defined by 

^■ = <^(0),Hess/(^(0))>. 

Note that ^(0) = ^ if V/(x) ^ 0. 

As in the proof of Theorem 13. 1\ we want to show that the solution 
to the Riccati equation 



(21) 



S t + S 2 t +Rt = 0, S = S 



is bounded below under the assumptions of the theorem. 
We compare the equation in (T2~Tj) with 



(22) 



S + S 2 + 




k\Vf\ 2 I 



0, So = S. 



By Theorem 13.41 and (1201) . we have S t > St- If we apply Theorem 
13. 5\ the solution of the initial value problem in (122 j) is given by 



5 f = r 1 (t)(r a (*))- 1 



where A 



ri(t) 



l*liv/|, 







cosh(At)/ 



1 

cos(At)/ 



5 



5 





Asinh(At)/ 


Asin(At)/ 



fc < 
k = 
> 
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and 



r a (f) 





sinh(At) 



ts + i 



t 





A 







sin(At) 
A 



Therefore, our assumption 

-1 



S > I 



I 







s 



s 



1 

cosh(At)/ 

1 

cos(At)/ 







-Acoth(A)/ 





-Acot(A)/ 



k < 
k = 
fc > 0, 



A; < 
k = 
A; > 0. 



implies that S'j is bounded for all t in [0, 1]. 

The rest of the proof follows as in the proof of Theorem 13 . 11 



□ 
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